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Abstract 



A two-site double exchange model with a single polaron is studied using a perturbation 
expansion based on the modified Lang-Firsov transformation. The antiferromagnetic 
to ferromagnetic transition and the crossover from small to large polaron are inves- 
tigated for different values of the antiferromagnetic interaction (J) between the core 
spins and the hopping (t) of the itinerant electron. Effect of the external magnetic 
field on the small to large polaron crossover and on the polaronic kinetic energy are 
studied. When the magnetic transition and the small to large polaron crossover co- 
incide for some suitable range of J/t, the magnetic field has very pronounced effect 
on the dynamics of polarons. 
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1. Introduction 



The origin of ferromagnetism in manganese perovskites Lai-xX^MnO^ {X = Ba, 
Sr, Ca etc.) is the double exchange mechanism [|T],3. The discovery of anomalous 
magnetotransport phenomena |^ in these compounds has stimulated intensive studies 
in its magnetic as well as electrical properties. However, the simple double-exchange 
interaction alone is not sufficient to explain the experimental results P]. Depending 
on the doping, temperature and the radius of the dopant ion these oxides show various 
phases 0. A complete understanding of the properties of Mn perovskites is still far 
from clear. 

One of the keys to understand unusual physical properties is to find the role of 
coupling between the carriers and the underlying lattice. Several theoretical models 
have been proposed based on lattice-carrier coupling p|-p!0[]. Many recent experiments 



||Tl| indicate that the electron-phonon (e-ph) interaction shapes its properties very 
crucially. Moreover, small to large polaron crossover is reported by many experi- 
mental groups |T^-|T5[. There are models P,|7|,pi][] which incorporate double exchange 



interaction in a polaronic model. Min and co-workers have studied the role of 
e-ph interaction in systems where double exchange interaction is present. The com- 
bined model of spin double exchange and lattice polaron is used to investigate the 
effect of small to large polaron crossover on the magnetic and transport properties 
under the mean-field approximation scheme. It concludes that the effect of polaron 
narrowing on the colossal magnetoresistance (OMR) is more pronounced in OMR 
manganites with low magnetic transition temperature (Tc) than in the high Tc man- 
ganites. Within the mean field theory the magnetic transition, metal-insulator 
transition and a large drop in magnetoresistance occur at the same temperature. 

In an attempt to have a clearer view from a nearly exact calculation, we include 
double exchange interaction in a two site one polaron model and follow a perturbation 
expansion |]l6l based on a modified Lang-Firsov (MLF) phonon basis where the lattice 



distortions produced by the electron are treated as variational parameters |17-19 • 



This method shows good convergence as well as nearly exact results for almost 



the entire range of e-ph coupling for t/uj^ < 1, where ujq is the phonon frequency. 
We investigate the ferromagnetic (FM) to antiferromagnetic (AFM) transition, the 
crossover from large to small polaron, the behavior of the effective hopping of the 
itinerant electron and the kinetic energy as a function of the e-ph coupling strength 
for the ground state of the system. Effect of the magnetic field on the large to small 
polaron crossover and on the polaronic kinetic energy are also studied. 

The paper is organized as follows. In Sec. 2 we define the model hamiltonian 
describing different interactions and calculate different physical quantities which in- 
dicate the behaviour of small-large polaron crossover and AFM-FM transition as a 
function of e-ph coupling. In Sec. 3 we present the results obtained in our calculation 
and discussions. Sec. 4 contains the concluding remarks. 

2. Formalism 



The Holstein model including double exchange interaction for the two-site 



single-polaron system is described by the Hamiltonian 

Q 

H = ^eui^ - ^tcos(-)(cl^C2a + 4aCi^) 



t <ij> 

J 



where z =1 or 2, denotes the site. Cia (cj^.) is the annihilation (creation) operator 
for the itinerant electron with spin a at site i, rii^ [=c\^Cia) is the corresponding 
number operator and g denotes the on-site e-ph coupling strength, hi and h\ are 
the annihilation and creation operators respectively for the phonons corresponding to 
interatomic vibrations at site i, ujq is the corresponding phonon frequency and e is the 
bare-site energy. Si and Sj are the core-spins at the site i and j respectively, 9 is the 
angle between the core-spins Si and Sj. The transfer hopping integral {t) is modified 



through the relative angle 9 a.st 003(2) because of the strong Hund's coupling between 
the spins of the core eletrons and itinerant electron J is the superexchange 
antiferromagnetic interaction between the neighbouring core-spins. Since we will 
restrict ourselves to the single-electron case we will not consider the electron spin 
indices. 

If one compares our simple model with the manganite (e.g. Lai^^C a^MnO^) sys- 
tem the core spin and the itinerant electron may be identified with the t^g (localized) 
electrons and the (mobile) electron of Mn^+ of the manganite system respectively. 
It may be mentioned that for manganites the Jahn-Teller(JT) coupling is important 
whereas the Holstein phonon mode, in general, should correspond to the breathing 
mode in manganites |Q. To include JT coupling one should consider two eg orbitals 
(say a and (3) and their coupling to the JT phonon mode. For the single electron 
case, as considered here, if one neglects the interorbital hopping {taj3 = 0) then the 
Hamiltonian can be separated into two parts each corresponds to 1-orbital case and 
as a consequence the 2-orbital problem is reduced to an effective single orbital prob- 
lem with an effective Holstein type interaction as in Hamiltonian (1). For the two 
electron case the above simplification cannot be done and a study of two-site two- 
electron problem considering two orbitals in the context of manganites is in progress. 
We have considered here one itinerant electron in a two-site 1-orbital system. So, 
the density of itinerant electron < n > is equal to 0.5 which, in principle, corresponds 
to X = 0.5 for manganite system. It may be mentioned that the studies of 1-orbital 
and 2-orbital Kondo model in the context of manganites using Monte Carlo 
techniques showed that the results for 2-orbital model at n = 1 (x = 0) are similar 
to those for single orbital model at x = 0.5. 

By introducing new phonon operators a = (&i + 62)772 and rf= (61 -62)772, the 
Hamiltonian (H) is separated into two parts : one corresponding to the in-phase mode 
which does not couple with the electronic degrees of freedom and the other involving 
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out-of-phase mode (Hd) which represents an effective e-ph system and cannot be 
solved analytically 

The MLF transformation with variable phonon basis is used so that a convergent 
perturbation expansion can be obtained. The transformed Hamiltonian is 

= uoSd + ^ epUi — t cos(-) [c\c2 exp(2A((i^ — d)) 

i 

+ 4ci ex.^{-2\{d^ -d))\+ujQ{g+-\){ni-n2){d + d^) + JS'^cosO (2) 

where R = X{ni — n2){d''—d), A is a variational parameter related to the displacement 
of the d oscillator, g+ = g/ and ep = e — LUQ{2g^ — A) A. 

For the perturbative expansion following Ref. |]16[ the basis set is chosen as 
|±,A^) = -^(ci ± C2) |0)e|iV), where |+) and |— ) are the bonding and antibond- 
ing electronic states and |A^) denotes the A^th excited oscillator state within the MLF 
phonon basis. The diagonal part of the Hamiltonian Ha in the chosen basis is treated 
as the unperturbed Hamiltonian (Hq) and the remaining part of the Hamiltonian 
Hi = Hd — Hq, as a perturbation. 

The unperturbed energy of the state |±, N) is given by 



E^^^^ = {N,±\Ho\±,N) 



NUo + tpTteff 



u=o ^■ 



+ JS^ COS 9 (3) 



where tg// = i cos | exp (— 2A^) and the general off-diagonal matrix elements of Hi 
between the two states |±, A^) and |±, M) may be calculated for (A^ — M) > as in 



Ref. \16 



Within the chosen basis, the unperturbed ground state is the |+)|0) state with 
the unperturbed energy, Eq^^ = ep — t^ff + JS"^ cos 6. 

The first order correction to the ground state wave function is obtained as, 

I ,(l)^ _ - A) - 2 Ate//] |_ .\ V- teffi2X)^ 



(E^-EZ) N=r,A..^N\{Er -E, 



e,NJ 



where e = + or - for even and odd N respectively. 

The first order correction to the energy {Eq'^) is zero since Hi has no diagonal 
matrix element in the chosen basis. The second order correction to the ground state 
energy is given by 

(2) _ ^ I ^^ + UJo{g+ - X) dN,l\ 

2^ ,p(0) p(0) X 

N=l,3.. l-C/Q - -C/_,Arj 

| -W/(2A)^ |2 

Higher order corrections to the ground state wavefunction and energy are obtained 



following Ref. |[T^. For the study of the effect of the magnetic field (h) we include a 



term -gJ2i fJ'shSi cos | to the Hamiltonian in Eq. (1) (where g is the Lande g factor) 
and as a result a term — 2/ie//^cos| is added to -Eq^'* where /ig// (= qS^b) is the 
local moment of the core spins. In this paper we express the magnetic field {h) in a 
unit of jieff = 1- 

Now a proper choice of A is to be made so that the perturbative expansion becomes 



convergent. Our previous work ||T6[ has shown that the A, obtained by minimizing 
the unperturbed ground state energy, gives satisfactory convergence to the perturba- 
tion series for t/ujQ < 1 and the convergence becomes very rapid and excellent with 
decreasing value of t/ujQ. Here also we will follow the similar procedure. Minimizing 
the unperturbed ground state energy E^^^^ with respect to A we obtain 

A = (6) 

Minimization of the unperturbed ground state energy with respect to 6 gives an 
approximate value (^mlf) of 9 

cos = [ h (loi' nonzero solution of 6'mlf) 

However, exact 9 should be evaluated from the minimization of the exact ground 
state energy. For each value of g+ we calculate the energy upto the sixth order in 



perturbation |T6| and find out for wliich value of 6 the energy (including perturbation 
corrections) is minimum. 

We have also calculated = —Exin =< V'gI ^cos(|) [c|c2 exp {2X{S — d)) 
+C2C1 exp(—2X{S — d))] \ipQ > where ipc is the ground state wave-funcion which we 
have calculated upto the fifth order corrections in perturbation, tfj^ describes the 
kinetic energy of the system and it reduces to the effective hopping (te//) for both 



the weak and strong coupling limits [22 . 

The static correlation functions {niUi)o and {niU2)o, where Ui and U2 are the 
lattice deformations at sites 1 and 2 respectively, produced by an electron at site 1, 
are the standard measure of polaronic character. The correlation functions for this 



two-site system may be written as 

(niMi)o = ^ -{g+ + A) + 
(niM2)o = ^ -{9+ - A) 

where Aq = (^/'g|^i('^ + c^^)|^g) 



Ng 
Ng 



(7) 



where Ng is the normalization factor to the wavefunction. The variation of the 
physical quantity Xcorr/g+ = — < "^1(^1 — U2) >o / 9+ with the e-ph coupling 
strength manifests the nature of large to small polaron crossover. Its value becomes 
1 in the extreme small polaronic limit where the lattice distortion is very local and 
decreases significantly in the small to large polaron crossover region as g+ decreases. 



3. Results and 

In the combined model of polaron and 
of variation of 9, tfP. and Xcorr/g+ with (7+ 



Discussions 

double exchange interaction, the nature 
reveal different transitions and crossover 



regions. We focus our attention to the parameter space where t > J as it corresponds 
to real systems [^,23|. We present the results for t = 1.0 and JS"^ = 0.05,0.15, 0.25 
and 0.5 (in a scale of ujq = 1.0). The bare-site energy e is taken as zero. 

Fig. 1 shows the change of relative orientation (6) of two core spins as a function 
of e-ph coupling strength (g^) for the ground state. For small values of two core 
spins are either aligned parallel or canted depending on the strength of t/JS'^. With 
increasing g^ the angle {6) increases and finally a transition from the FM or canted 
AFM state to the AFM state occurs. Nature of this transition depends on the ratio 
of t to JS'^. For t=l.O, JS^=0.05 the FM to AFM transition is very sharp. With 
increasing value of JS"^ the transition occurs at a lower value of g^ and it becomes 
broader. For higher values of JS"^ (=0.25, 0.5), the crossover from the canted state 
to the AFM state is very smooth. It may be mentioned that Yunoki and Moreo |§ 
studied the 1-orbital Kondo model with AF interaction (J) between the core spins 
(in absence of e-ph interaction) in the context of manganites and obtained FM state 
ai n = 0.5 for low values of J. With increasing J the ground state at quarter filling 
(n = 0.5) appears to be a spin spiral state and ultimately an AF state. For the two 
orbital model in presence of e-ph coupling a similar FM state is observed for low 
values of e-ph coupling for n = 1 and an AF state for large values of e-ph coupling. 

With increasing g+ a large to small polaron crossover occurs in the system in 
addition to the magnetic transition. This is clearly seen from Figs. 2 and 3 where 
the variations of t^/, ^corr/g+ and 6 with g^ are shown. The large to small polaron 
crossover is identified by the sharp fall in and rise in Xcorr- For t = 1.0 and 
JS"^ = 0.15 the FM to AFM transition and the large to small polaron crossover 
occur simultaneously (Fig. 2) around g+ = 1.2. However, for JS^=0.05 the magnetic 
transition occurs at a higher values of g^ than that coresponding to large to small 
polaron crossover (Fig. 3). In this situation the FM state has two regions: one with 
appreciable value of and the other with much reduced value of tfj^ (Fig. 3). 
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In the former case, charge carrier is the delocahzed large polaron with appreciable 
hopping and this regime is expected to be metallic in the thermodynamic limit. In the 
latter region the hopping would be heavily suppressed by the small polaron formation 
and it may result in an insulating FM state in some cases. Thus we find that there 
is a possibility of FM metal-FM insulator- AFM insulator transition for low values of 
JS'^/t. It is to be noted that if one follows the usual MLF method (zeroth order of 
perturbation within our approach) the FM-AFM transition coincides with large to 
small polaron crossover for any value of JS'^/t. This is demonstrated in Fig. 3 for 
JS"^ /t=0.05, where we find the coincidence of the FM-AFM transion with the large- 
small polaron crossover within the usual MLF approach whereas the nearly exact 
calculations yield different result. The nearly exact results, obtained from the MLF 
perturbation method, makes the large-small polaron crossover smooth in favour of 
the conclusion of Lowen . 

In Fig. 4 we have shown the effect of the external magnetic field on 6, Xcorr/g+ and 
tjf for JS'^ = 0.15. The magnetic field favours the FM state, hence the FM-AFM 
transition and associated large-small polaron crossover (for JS^/t = 0.15) take place 
at a higher value of with increasing field. From Fig. 4 it is clear that the magnetic 
field reduces the value of the local distortion {Xcorr/g+) and enhances significantly 
in the large-small polaron crossover region provided the crossover is associated with 
a FM-AFM transition. 

For JS'^/t = 0.05 the system remains in the FM state in the large-small polaron 
crossover region, consequently the magnetic field has no effect on in this crossover 
region. However, the field has an effect on 9, hence on t^fj at the FM-AFM transition. 

In Fig. 5 we plot the change in due to the magnetic field as a function 
of This quantity may be related to the magnetoresistance for a system in the 
thermodynamic limit. For polarons the conductiviy at low temperatures is dominated 
by the tunneling mobility which is proportional to t^fj within the zeroth order of 
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perturbation |10]. In general the is a measure of delocalization of the electron. A 
reduction of will cause a reduction in mobility. In Fig. 5 the change in if^^ due to 
the field is very sharp and prominent at the FM-AFM transtion when it is associated 
with the large-small polaron crossover for JS'^/t = 0.15, while it shows a broad 
smaller peak for JS'^/t = 0.05 where the FM-AFM transition is not associated with 
the polaron crossover. For the former case the magnetic field would have pronounced 
effect on the transport properties. 



4. Conclusions 



From our studies on the two-site double exchange model with a single polaron as 
a function of e-ph coupling strength we conclude that the nature of the FM-AFM 
transition depends on the relative values of J and t. The transition is sharper for 
smaller values of J/t. For high values of J/t a canted state is stable instead of a 
FM state for weak e-ph coupling ((yf+) and the crossover from the canted state to the 
AFM state is very smooth with increasing g^. For suitable values of J/t (=0.15) 
the FM-AFM transition coincides with the large-small polaron crossover. For this 
case the external magnetic field has very prominent effect on the polaronic local 
distortion and kinetic energy in the transition region (as in Figs. 4(b) and 4(c) 
respectively). Within the MLF method (zeroth order of perturbation) the polaron 
crossover is always associated with FM-AFM transition for any value of J/t whereas 
our results based on convergent perturbation expansion yield that the coincidence 
of magnetic transition and polaron crossover will depend on the value of t and J. 
For low values of J/t a crossover from the FM large polaronic state to the FM small 
polaronic state with reduced hopping occurs and then a transition to the AFM state 
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takes place at a higher value of 
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Figure captions : 

FIG. 1. Variation of the relative angle (9) between the core-spins, with g+ for 
t = 1.0 and JS'^ = 0.05, 0.15, 0.25 and 0.5 (in units of cuo=l). 

FIG. 2. Variations of 9, tf^f and \corr/g+ with g+ for t = 1.0, JS^ = 0.15 and 

FIG. 3. Variations of 9 (with and without perturbation), tff^, tgff, Xcorr/9+ ^-nd 
X/g+ with g+ for t = 1.0, JS"^ = 0.05 and h = 0. Note that ^(MLF), te/f and \/g+ 
are obtained within the MLF method without considering perturbation corrections. 

FIG. 4. Variation of (a) 9, (b) Xcorr/g+ and (c) t^f^ with gf+ for different values of 
the magnetic field h =0, 0.02 and 0.04 (in units of /jL^ff^l) for t = 1.0 and JS^ = 0.15. 

FIG. 5. Variation of the change in igf due to the magnetic field (l!llM^!liM) 
as a function of g+ for t = 1.0, JS^ = 0.05 and 0.15 and /i=0.04. 
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